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Abstract 

We consider the relation between the Q deformed N = 2 SQCD with the single 
deformation parameter e and integrable models of the BCS-like superconductivity. It 
is argued that the vortex string worldsheet theory is related to the Russian Doll(RD) 
model of the truncated BCS superconductivity. We argue that the Q deformed gauge 
theory manifests the interesting cyclic RG behavior with the period of the RG cycle 
proportional to e _1 . The deformed gauge theory can develop several non-perturbative 
scales. We conjecture on the monopole bound state interpretation of the Efimov 
tower. 



1 Introduction 



A confinement of electric degrees of freedom implies a condensation of some magnetic 
degrees of freedom. They certainly emerge from D-branes however the precise pattern 
of the D-brane condensation responsible for confinement in the asymptotically free 
gauge theories is still unclear. One could imagine the condensation of the monopoles, 
strings, domain walls or complicated webs designed from these ingredients. Therefore 
it is useful to look at the situations where the gap generation could be understood 
from the first principles. The simplified models of superconductivity provide the 
proper background to analyze the different aspects of these issues. 

The well-known Seiberg-Witten solution [T] yields the pattern for the monopole 
condensation in the perturbed N = 2 SUSY YM theories with the different matter 
content. It was argued that the low-energy effective action and the spectrum of the 
stable BPS states are governed by the classical finite-dimensional integrable system 
[2]. In particular the theory which we shall mainly focus on, namely Nf = 2N C SQCD, 
is governed by the twisted XXX spin chain [3]. The vacua in the deformed theory 
correspond to the classical equilibrium states in the corresponding integrable system. 
The relation between the finite dimensional integrable systems and supersymmetric 
YM theories has been generalized to the quantum case [I]. The parameter of the 
Q deformation plays the role of the Planck constant and the Bethe anzatz equation 
in the finite dimensional integrable system corresponds to the extremization of the 
superpotential in the deformed gauge theory. 

On the other hand the finite-dimensional integrable systems in some cases pro- 
vide the exact solution to the truncated BCS - like models of superconductivity. The 
generic arguments explaining the relation between the integrable systems and the 
superconducting models go as follows. Consider the fermionic system with a kind of 
attractive four-fermion interaction. It can be derived upon the integration over the 
"phonon"-like degrees of freedom. The system develops the fermionic BCS conden- 
sate of the Cooper pairs which is generically inhomogeneous. It is assumed that the 
single fermion or Cooper pair propagating on the top of the condensate should not 
ruin it. This is a kind of consistency condition for the combined system " inhomoge- 
neous condensate + excitation" . It is this consistency condition that provides some 
integrable system behind the scene. The space-time dependence of the condensate 
is governed by the integrable equation itself while the equation for the one-particle 
fermion wave function in the condensate background is interpreted as the Lax equa- 
tion for the corresponding dynamical system. The condensate can be represented 
via the continuous gap function or some lattice of Abrikosov strings for the type II 
superconductivity. In the first case one has the continuous integrable system of KdV 
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type while in the second case the discrete systems of Toda type are involved. 

The analogy with the Peierls superconductivity has been developed in j5] for the 
pure perturbed N = 2 SYM theory whose low-energy sector is described by the 
periodic Toda chain [2]. The Pieirls model involves the fluctuating one- dimensional 
crystal and fermions propagating in such background [Bj. The Toda system gov- 
erns the interaction between phonons induced by fermions while the Lax equation 
describes the propagation of the fermion in the background of the emerging conden- 
sate. The Seiberg-Witten curve is identified with the dispersion law for the fermions 
and the nonperturbative Aq C d scale in SYM corresponds to the superconducting 
gap in the Peierls model. 

The Toda system emerges as the particular limit of the XXX spin chain when the 
Casimirs of the representations at each site of the chain tend to infinity. Hence we 
could expect some superconducting model behind the SQCD as well generalizing the 
arguments from (5j. It turns out that such model exists indeed and we shall describe 
it in the paper. The Richardson model [7] and the RD model [8] are proper gener- 
alizations for pure SQCD and e-deformed SQCD correspondingly. The Richardson 
model can be also considered as the XY model in the inhomogeneous external field. 

On the other hand the Richardson model links up with the generalized Gaudin 
system [9] and the RD model with the twisted inhomogeneous XXX spin chain 
correspondingly [10]. It turns out that the relation between the quantum integrable 
systems and the vacua of the supersymmetric gauge theories [I] is quite convenient 
to clarify the underlying physical degrees of freedom. The Bethe anzatz equation 
describing the spectrum of the RD model coincides with the equation describing the 
ground state in the worldsheet theory of several nonabelian strings. These semilocal 
nonabelian strings with rich worldsheet theory [Til EG] • are the classical solution to 
the equations of motion in the SQCD . 

To some extend the integrability plays the role of consistency condition in more 
general setting. Recall that the most familiar example of such phenomenon is pro- 
vided by the strings propagating in the external metric or the gauge field. It is 
well-known that the condition of the vanishing of Fradkin-Tseytlin beta-function 
of the string worldsheet theory corresponds to the classical equations of motion for 
background fields. Speaking differently this is the consistency condition for the prop- 
agation of the string in its own condensate. 

In the problem under consideration the situation is more subtle since the world- 
sheet theory of the nonabelian string is not conformal. Hence the consistency con- 
dition is more involved and the quantum properties of the bulk fields have to be 
matched in this case. The interplay between the D = 4 bulk theory and the D = 2 
worldsheet theory on the nonabelian string has been examined in [121 Q3] where it 
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was argued that the /3-functions and the spectrum of the BPS states match. In 
this paper we shall argue that the additional ingredient can be added to the generic 
matching condition that is the system of superconducting fermions. The correspond- 
ing Richardson of RD models are conjectured to describe the magnetic degrees of 
freedom and the consistency condition implies that their dispersion law fits with the 
spectral curve of the proper integrable system. These integrable spin chain models 
correspond to the dynamics of the nonabelian strings per se. 

Due to the mapping of e -deformed SQCD into the RD model it is possible to 
gain interesting observations from the known results concerning the superconducting 
model. The most unusual feature of RD model is the cyclic RG behavior yielding 
the multiple condensates instead of the single one [8]. The period of the RG cycle 
is determined by the RG invariant and the spectrum of the model is reshuffled upon 
each cycle. There are several systems enjoying cyclic RG flows both in the quantum 
mechanics [15] and in the field theory with the several couplings [16]. In particular 
such behavior has been found in the sin-Gordon model analytically continued in the 
coupling constant [16] . In that example the cyclic RG behavior amounts to the set of 
the unusual resonances with the Regge-type stringy spectrum. In general the origin 
of the cyclic RG behavior is some resonance-like behavior in two-body system which 
amounts into the hierarchy of the Efimov-like states. The most surprising aspect of 
the cyclic RG is its sensitivity to the UV scale of the theory under consideration. 
Note that is was shown recently that the cyclic RG flow is consistent with the c- and 
a-theorems [PT] . 

The RG step in RD model corresponds to the change of the XXX spin chain 
length. Being translated to the conformal gauge theory side it corresponds to the 
decoupling of two flavors sending the corresponding mass to infinity with the simul- 
taneous change iV — > (N — 1). The theory remains conformal with the different rank 
of the gauge group. However in the deformed theory we have additional dimension- 
less parameter which is the ratio — hence decoupling of the heavy flavor is described 
via two-coupling RG flow. Using the known results concerning RD model we shall 
argue that decoupling of a heavy flavor occurs in a cyclic manner and the period of 
a cycle is proportional to the inverse of the deformation parameter therefore in the 
undeformed case the period is infinite. 

The remarkable AGT relation [19] makes manifest the connection between the 
SU(2) Nekrasov partition function of the deformed SYM theory and the conformal 
blocks in the Liouville theory [19] . The interplay between the four-dimensional and 
two-dimensional theories reflects the complimentary viewpoints of observers at the 
M5 brane worldvolume where the SYM lives on. For the higher rank group the 
Liouville theory gets substituted by the 2d Toda theory [20] ■ The nonabelian strings 
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with large tension in the gauge theory are identified as the surface operators which 
on the other hand are represented as the proper vertex operators in the Liouville- 
Toda theory [21] . In particular the D2 branes representing the surface operators 
in Liouville-Toda theory yield the degrees of freedom for the Toda-Calogero type 
integrable models. 

Hence we could look at the possible interpretation of the cyclic RG flows for 
the conformal blocks in the Liouville/Toda theory and the possible interpretation of 
the superconducting models on the gravity side of the AGT correspondence. The 
key point is the interpretation of the fermions and their Cooper pairs from the first 
principles. We shall make some conjectures however the complete analysis is still to 
be done. 

The paper is organized as follows. In Section 2 we briefly review the Peierls 
model while in Section 3 the Richardson model and the RD model are considered. In 
Section 4 we comment on the generic properties of the cyclic RG flows. In Section 5 
the mapping to the vortex nonabelian states in the SQCD is considered. In Section 
6 we conjecture on the interpretation of the fermions in the superconducting model 
as the monopoles in the Higgs branch. The main observations of the paper and the 
list of the open questions can be found in the last Section. 

2 The Peierls model 

In this section we briefly review the main facts concerning the Peierls model relevant 
for the further consideration. Originally it was formulated to describe the selfcon- 
sistent behavior of Id fermions interacting with the fluctuating lattice and has been 
applied to the description of the Id superconductivity [6]. The Coulomb interaction 
between the fermions is neglected, the fermions are assumed to be in the external 
field determined by the lattice degrees of freedom while the lattice dynamics itself 
gets modified by the fermions. In what follows we will mention the integrable struc- 
ture both in the continuum [22] and discrete Peierls models [6] . Hamiltonian density 
for the simplest continuum model looks as follows 

H con = ^ + a 3 d x ^ + ^+(cx_A(a;)* - a + A(x))V + A(x) 2 (2.1) 

where A (a;) represents the lattice potential. The saddle point solution for A(x) 
provides the four-fermion interaction Hamiltonian induced by the phonon exchange. 
To some extend the lattice potential A(x) can be interpreted as the inhomogeneous 
fermionic condensate. 
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For the discrete version one has [6] 

H dis = ^(*+V„* n + V+C n V n+1 + *+C„_i* n _i) + ^ Ki^i (2-2) 
n i 

c n = exp(x n+ i - x n ), I N = ^ lnc n, hN = + v l) ( 2 -3) 

n n 

where x n ,v n are the lattice coordinate and momentum while I n are identified as the 
Toda chain Hamiltonians. 

To identify the ground state of the model we minimize the Hamiltonian with re- 
spect to the fermionic and lattice variables. The variation over the fermionic variables 
yields the Lax equation for the Toda chain 

C n ^ n+i + C„_i* n _i + V n ^! n = E^ n (2.4) 

Variation over the lattice degrees of freedom gives rise to the system of the finite 
number of algebraic equations describing the Riemann surface. In the Peierls model 
we assume the periodicity of the fermionic wave function on the lattice 

* n+N (E) = e iN ^ n (E). (2.5) 

The solution is described in terms of the hyperelliptic Riemann surface 

y + y- 1 = P N (.E) (2.6) 

where P^(E) is the polynomial. This Riemann surface plays the double role. First, 
its Jacobian is identified as the complex Liouville torus for the periodic Toda chain 
describing the dynamics of phonons. On the other hand it is nothing but the Fermi 
surface for the fermions. The moduli of this surface are fixed by the minimization 
equations in terms of the physical parameters like the fermionic density p = 
The key feature of the solution is the appearance of the fermionic mass gap which is 
analogue of the Aqcd- Fermionic wave function is uniquely defined on this surface 
and the number of its zeros coincides with the genus of the curve. It was proved 
that only one gap vacuum configuration is stable, therefore the generic polynomial 
degenerates. The very similar one gap vacuum configuration emerges in the Gross- 
Neveu model [T8] . 

The spectrum of the excitations above the ground state involves fermionic and 
bosonic quasiparticles. The phonon and the charge density wave represent the 
bosonic excitations j22l [23] . The fermionic excitations strongly depend on the 
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fermionic density p. At large p the polyaron type state is localized at the lattice 
configuration 

AO) = const - ——^ — r (2.7) 

In the opposite limit at small density one has the delocalized fermionic state and the 
gap potential 

A(a?) = const + xcos(2y/xx + (p) (2.8) 

where x is some constant. Exact solution provides the temperature dependence of 
the mass gap. It was shown [22] that the fermion mass gap gets renormalized and 
disappears at some critical value T c . Being translated into the form of the dispersion 
law it tells us that the Riemann surface degenerates to a sphere above the phase 
transition point. 

Let us compare the data governing the vacuum structure of the pure N = 2 SUSY 
YM theory and the one from the Peierls model. Low energy effective action in SYM 
theory is fixed by the Riemann surface and holomorphic differential defined on it [I] . 
Prepotential T can be derived from the relations 

an = a; = J) A an = (f A A = dS = Edp (2.9) 

J Ai 1 J B . 

where and Bi are the cycles on the Seiberg-Witten spectral curve. The number 
of cites in the Peierls model corresponds to the rank of the gauge group while the 
total length can be identified with the coupling constant r o = ~? + taken at the 
UV scale. The density of the fermions turns out to be related to the renormalized 
coupling constant in the gauge theory [5]. The equations selecting the ground state 
of the Peierls model correspond to the selection of the point at the Coulomb branch 
of the moduli space. 

It what follows we shall generalize the relation between the low-energy sector of 
SYM theory and models of superconductivity for the theories with matter in the 
fundamental representation. 



3 Truncated models of BCS superconductivity 

3.1 Richardson model versus Gaudin with irregular singu- 
larities 

Let us recall the truncated BCS-like Richardson model of superconductivity [7] with 
some number of doubly degenerated fermionic levels with the energies e J(J . . . j = 
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1 . . . N. It describes the system of a fixed number of the Cooper pairs. It is assumed 
that several energy levels are populated by Cooper pairs while levels with the single 
fermions are blocked. The Hamiltonian reads as 

N 

H B cs= ^ e ja cf r c ja -G^2cj + cj_c k -c jk+ (3.10) 

j,a=± jk 

where are the fermion operators and G is the coupling constant providing the 
attraction leading to the formation of the Cooper pairs. It terms of the hard-core 
boson operators it reads as 

Hbcs = ( Ai h J ~ G J2 h P k ( 3 - U ) 

j jk 

where 

[bjb k ] = 5 jk (2Nj - 1), 6 i = c i _c,-+, N 3 = bp, (3.12) 
The eigenfunctions of the Hamiltonian can be written as 

M N 

\M>=l[B t (E l )\vac>, = — 6+ (3.13) 

i 3 ■> 

provided the Bethe anzatz equations are fulfilled 




(3.14) 



The energy of the corresponding states reads as 

E{M) = Y,Ei (3.15) 

i 

It was shown in [9j that the Richardson model is exactly solvable and closely 
related to the particular generalization of the Gaudin model known as the model 
with irregular singularities [26J. To describe this relation it is convenient to introduce 
the so-called pseudospin si (2) algebra in terms of the creation- annihilation operators 
for the Cooper pairs 

r = b t + = b + t° = N -1/2 (3.16) 
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The Richardson Hamiltonian commutes with the set of operators Rj 



t>it>j 



R i = -t° i -2G'52-!¥i- (3.17) 



which are identified as the Gaudin Hamiltonians. 



[Hbcsi Rj] = [R i ,R j }=0 (3.18) 

Moreover the Richardson Hamiltonian itself can be expressed in terms of the opera- 
tors Ri as 

hbcs = e ^ + G E R ^ 2 + const ( 3 - 19 ) 

3 

The number of the fermionic levels N coincides with the number of sites in the 
Gaudin model and the coupling constant in the Richardson Hamiltonian corresponds 
to the "twisted boundary condition" in the Gaudin model. The Bethe anzatz equa- 
tions for the Richardson model ( 13.141) exactly coincides with the ones for the gen- 
eralized Gaudin model. Note that the Gaudin model is the Hitchin system on the 
marked sphere which plays the role of the Fermi surface for the fermions very much in 
the same manner as we have seen in the Peierls model. It was argued in [8] that the 
Bethe roots corresponds to the excited Cooper pairs that is natural to think about 
the solution to the Baxter equation as the wave function of the condensate. In terms 
of the conformal field theory Cooper pairs correspond to the screening operators [25J. 

For the nontrivial degeneracies of the energy levels dj the BA equations read as 




(3.20) 



3.2 Russian Doll model of superconductivity and twisted 
XXX spin chains 

The important generalization of the Richardson model describing the reduced BCS 
superconductivity as well is the so-called RD model [8]. It involves the additional 
dimensionless parameter a and the RD Hamiltonian reads as 

N 

H RD = 2 ^(e, - G)N t - G ^{eH+ty + e~ ia b+b k ) (3.21) 

i j<k 
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with two dimensionful parameters G, rj and G = \/G 2 + rf. In terms of these vari- 
ables the dimensionless parameter a has the following form 

a = arctanh(-) (3.22) 

Gr 

It is also useful to consider two dimensionless parameters g, 9 defined as G = gd and 
Tj = 9d where d is the level spacing. The RD model reduces to the Richardson model 
in the limit 77 — >■ 0. 

The RD model turns out to be integrable as well. Now instead of the Gaudin 
model the proper counterpart is the generic quantum twisted XXX spin chain [10] 
The transfer matrix of such spin chain model t(u) commutes with the Hre> which 
itself can be expressed in terms of the spin chain Hamiltonians. 

The equation defining the spectrum of the RD model reads as 



exp(-2ia) TT § ~ £fc ~ ^ 2 = ff e i- e *-*7 ( 3 .23) 
Pl >lL Ej -e k +iri/2 1 1 e j - e k + irj 



and coincides with the BA equations for the spin chain. It reduces to the BA equation 
of the Richardson model (13.141) in the limit 77 — ?• 0. 

The key feature of the RD model is the multiple solutions to the gap equation. 
The gaps are parameterized as follows 

u Tin , ,. 

n = ImM' t n = t + — n = 0, 1 . . . (3.24) 

where i is solution to the following equation 

9 7T 

tan(6t ) = - < t < 77 (3.25) 
9 

and bj = dN for equal level spacing. This behavior can be derived via the mean 
field approximation [21]. Different solutions to the gap equation yield the different 
superconducting states. In the limit 9 — > the gaps A n>0 — » and 

t = -, A = 2wei (3.26) 
9 

therefore the standard BCS expression for the gap is recovered. At the weak coupling 
limit the gaps behave as 

A n oc A e-ir (3.27) 

In terms of the solutions to the BA equations the multiple gaps correspond to the 
choices of the different branches of the logarithms. 
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If the degeneracy of the levels is d n then the RD model gets modified a little bit 
and is related to the higher spin XXX spin chain. The local spins Sj are determined 
by the corresponding higher pair degeneracy di of the i-th level 

Si = di/2 (3.28) 

and the corresponding BA equations read as 

A Ej -e k -ir}/2 + irjSj _ -^r 

Ej - e k + ir]/2 - ir]Si H e 3 - - e k + ir] 



When the local spin tends to infinity the twisted XXX chain degenerates into the 
periodic Toda chain hence the large degeneracy of the levels is necessary at the 
superconducting side. It is this limit which corresponds to the Peierls model discussed 
earlier. 



4 Cyclic RG flows and RD model 

Let us explain the key points concerning the cyclic RG phenomena. It was anticipated 
for a while however the first clear-cut example has been elaborated only in [15] for 
the finite-dimensional system. The method developed in [15] formulates the step in 
the RG flow in the finite-dimensional system as removing the single highest energy 
level with simultaneous renormalization of the couplings. This approach has a lot 
in common with the renormalization procedure in the matrix models considered in 
[28] . The different approach to RG in the finite-dimensional system concerns the 
introducing the UV cutoff at the small scales and looking at the cutoff dependence. 
This approach has been successfully applied to the rational Calogero model [27] and 
it was argued that conformal quantum mechanics manifests the cyclic RG which 
reflects a kind of " anomalous" violation of the conformal group down to the discrete 
one. 

In fact the important phenomena which has the cyclic RG origin has been dis- 
covered long time ago by Efimov in the three-body system . His key observation 
concerns the emergence of the specific bound states in the three-body system involv- 
ing two very different scales. It was argued that in such situation when there are the 
two-body bound state near threshold there is the tower of the bound states with the 
RD scaling behavior of energies 

E n+ i = e s E n (4.30) 
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The presence of so-called Efimov states with the RD scaling is the common feature of 
the models with cyclic RG phenomena. The spectrum is reorganized in the universal 
way during the single cycle and the total number of Efimov states is of order log( 1 - ILE -). 
The review on the cyclic RG in the finite-dimensional systems can be found in [29|. 

There are examples of the cyclic RG flows in the 2D field theories with the several 
couplings where the RG is formulated in the standard manner as the dependence on 
the log of the renormalization scale [16]. It was argued that generically there should 
be at least two couplings and usually one of them does not run at all. In the field 
theory one also has the RD scaling for the Efimov-like resonances which in some 
examples manifest the Regge like structure. Moreover it was argued that the S- 
matrix behaves universally under the cyclic RG flows. The total number of Efimov 
states scales in the same manner as in the quantum mechanical case. 

The RD model of truncated superconductivity enjoys the cyclic RG behavior [8]. 
The RG flows can be treated as the integrating out the highest fermionic level with 
appropriate scaling of the parameters using the procedure developed in [301 EE] ■ The 
RG equations read as 

9n-i = 9N + ^(g 2 N + 2 ), 6 N -i = 6 N (4.31) 

At large N limit the natural RG variable is identified with logN and the solution to 
the RG equation is 

g(s) = 6tan(6s + tan" 1 ^)) (4.32) 

9 

Hence the running coupling is cyclic 

g(s + \)=g(s), g(e~ x N)=g(N) (4.33) 



with the RG period 



A = ~ Q (4.34) 
and the total number of the independent gaps in the model is 

Q 

N cond cc -logN (4.35) 

7T 

The multiple gaps are the manifestations of the Efimov-like states. The sizes of the 
Cooper pairs in the n-th condensates also have the RD scaling. The cyclic RG can 
be derived even for the single Cooper pair. 

What is going on with the spectrum of the model during the period? It was 
shown in [21] that it gets reorganized. The RG flows possesses the discontinuities 
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from g = +00 to g — —00 when a new cycle gets started. At each jump the lowest 
condensate disappears from the spectrum 

A n+ i(# = +00) = A n (g = -00) (4.36) 

indicating that (N+l)-th state wave function plays the role of N-th state wave func- 
tion at the next cycle. 

The same behavior can be derived from the BA equation [21]. To identify the 
multiple gaps it is necessary to remind that the solutions to the BA equations are 
classified by the integers mj , i = 1 , . . . M parameterizing the branches of the loga- 
rithms. If one assumes that rrii = m for all Bethe roots then this quantum number 
gets shifted by one at each RG cycle and was identified with the integer parame- 
terizing the solution to the gap equations At the large N limit the BA equations of 
the RD model reduce to the BA equation of the Richardson-Gaudin model with the 
rescaled coupling 

G' 1 = r ] - 1 (a + TTm) (4.37) 

which can be treated as the shifted boundary condition in the generalized Gaudin 
model parameterized by the integer. Let us emphasize that the unusual cyclic RG 
behavior in due to the presence of two couplings in the RD model. 



5 Nonabelian string versus BA equation 

5.1 BA equations at the string worldvolume 

Let us consider the Q deformed N = 2 SQCD with SU(L) gauge group, L funda- 
mental hypermultiplets with masses rrif. and L antifundamental hypermultiplets with 
masses m a fi. In the NS limit with the single deformation parameter e the Coulumb 
branch of the vacuum manifold parameterized by the vev of the adjoint scalar gets 
reduced to the set of points since the theory develops the twisted superpotential [I] 
determined by Nekrasov partition function. The isolated vacua are determined by 
the equation 

which yields 

a = rn f -N e (5.39) 

where n = (rii . . . n^). The integers can be attributed to the selection of the branch 
of logarithm in the twisted superpotential. 
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The bulk D=4 theory develops the Higgs branch where the scalar components 
of the fundamental are condensed. The Higgs branch admits the stable nonabelian 
string solution with the rich worldsheet theory. The physics of the non-abelian strings 
in the nondeformed case is described in the reviews [11] . 

In the bulk theory in the NS limit the worldsheet theory on the nonabelian string 
involves the L fundamental chiral multiplets with twisted masses Mp t and L anti- 
fundamental multiplets with twisted masses Maf v The additional chiral multiplet 
in the adjoint representation gets non-vanishing mass e due to the the background 
graviphoton field. It can be integrated out amounting to the twisted superpotential 
for the scalars in the vector multiplet. The minimization of the worldsheet theory 
at N nonabelian strings yields the BA equations where Bethe roots \ correspond 
to the values of the scalars in the vector multiplet [H [32], EI]. The solutions to 
the BA equations are parameterized by the set of integers rli obeying the condition 

It turns out that on-shall values of the twisted superpotentials in the bulk theory 
and the worldvolume theory of N nonabelian strings coincide 

W AD (cn =mi- rue) - W AD (ai = m x - e) = W 2D {n L ) (5.40) 

upon the following identifications of parameters [32] 

M f = m f - 3/2e, M af = rh af - l/2e (5.41) 

The rank of the worldvolume theory iV which corresponds to the number of non- 
abelian strings is defined via the relation 

N + L = J2 n h ni = m-l (5.42) 

i 

The modular parameter in the worldsheet theory 

t 2 d = ir + — — (5.43) 

A pi 

where r is the FI parameter related to the D = 4 coupling constant |13| . The 9 term 
penetrates the worldsheet theory from the bulk one |43j. The modular parameters 
of the theories are related as 

r 2D = n D + ^(N + l) (5.44) 

In the D = 4 bulk theory the BA equations emerge in the saddle point calculation 
of the instanton partition function [31]. The number of the Bethe roots coincides 
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with the number of the nonabelian strings which is in the perfect agreement with 
the interpretation of the nonabelian strings as the excitations above the root of the 
Higgs branch in the bulk theory. 

The BA equations for the deformed SQCD exactly coincide with the BA equations 
involved into the solution to the RD model. The asymmetry parameter of the RD 
model 77 is identified with the deformation parameter of the Q background 

77 = e (5.45) 

The fermionic energies in the RD model are identified with the masses of fundamen- 
tals in D = 4 bulk theory or the twisted masses in the worldsheet theory 

Ei = Mi (5.46) 

The modular parameter in the worldsheet theory provides the twisted boundary 
conditions in the BA equations for the inhomogeneous XXX spin chain. 

Geometrically the Bethe root or vev of the scalars corresponds to the position of 
the D2 brane representing nonabelian string in the transverse coordinates. Note that 
from the spin chain viewpoint the Bethe roots are zeros of the polynomial solution 

Q(X) = l[(X - Aj) (5.47) 

to the Baxter equation for the twisted spin chain. 



5.2 Cyclic RG in the deformed gauge theory 

Let us turn to the important observation. Since the BA equation for the deformed 
SQCD coincides with the one for the twisted XXX spin chain we can use the results 
concerning its cyclic RG behavior. First remind once again that there are two di- 
mensionless couplings in the Q deformed SQCD required for cyclicity. One of them is 
the conventional complexified coupling in the bulk gauge theory which plays the role 
of the twist in the spin chain. The second dimensionless parameter in the simplest 
case can be identified with the ratio 

9 = 4- (5-48) 
dm 

where 5m is the difference between the masses of the fundamentals. We assume that 
the masses are almost equidistant. It is useful following j2l] to introduce the second 
dimensionful parameter G similar to the Aqcd scale 

G 

tanr = — (5.49) 
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This scale is evidently nonperturbative with respect to the gauge coupling. 

To describe the cyclic RG behavior in SQCD let us identify the step of the RG 
flow in terms of the finite-dimensional system. As we have mentioned the step in the 
RD model corresponds to the decoupling of the highest energy level. Since the energy 
level in RD model corresponds to the mass of fundamental the RG step corresponds 
to the decoupling of the heavy flavor N F — > N F — 2. Simultaneously the rank of the 
group is changed N c — > N c — 1 and the form of the BA equation holds the same. 

In the non-deformed case this procedure yields the RG flow without any cycles. 
On the other hand if we just decouple the heavy flavor without changing the rank of 
the group the nonperturbative scale emerges in the asymptotically free theory. 

A = M ''^ pi ^MZM ] (5 ' 50) 

where M reg is the UV scale and Pq is the coefficient of the /3-function. In the de- 
formed case the situation is more involved since the decoupling of the heavy flavor 
is described by two-coupling RG equations. As we have seen before the e parameter 
is not deformed however the modular parameter enjoys the cyclic RG solutions. In- 
deed the complexified gauge coupling is expressed as a function of the scales via (??) 
hence the RG flow of G yields the RG flow of r. Let us emphasize that we consider 
the conformal theories and the cyclicity involves the rank of the gauge group. The 
ciclicity of the RG flow breaks the conformal group down to the discrete subgroup 
selected by the deformation parameter. 

As we have described above the key feature is the emergence of the multiple 
Efimov-like scales in the problem. In the RD model these Efimov-like scales corre- 
spond to the multiple gaps with the Efimov scaling. In the deformed SQCD this 
scales emerge as the multiple nonperturbative A-like scales whose number is defined 
by the ratio 

Ncond oc — (5.51) 

At the weak coupling the scales behave as 

A n oc A e — (5.52) 

which are certainly nonperturbative in the parameter of the Q background. 

Our observations provide the qualitative picture behind the RG flows in the de- 
formed QCD and more detailed analysis is required. In particular we have considered 
the situation with the single massive parameter only. In the generic situation we have 
the set of the dimensionless parameters 

6 % = ^ (5.53) 
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hence the RG flow involves the multiple couplings. 

6 On the Cooper pair interpretation 

So far we have not discussed the interpretation of the fermions developing the super- 
conductivity. In this Section we make some conjecture concerning their identification 
and present the several evidences supporting it. Namely we shall conjecture that the 
relevant degrees of freedom which form the Cooper pairs are the monopoles in the 
Higgs branch localized at the nonabelian string. 

6.1 On the monopole interpretation from the knot homolo- 
gies 

First let us make a few comments concerning the another appearance of the BA 
equations for the generalized Gaudin model. Attempting to get the field theory 
realization of the knot homologies the counting of the solutions to the BPS equations 
in the proper gauge theory was considered in [33]. The problem of counting of 
the BPS equations was reformulated in terms of the counting of the solutions to 
the BA equations in the particular integrable system. The mapping of the gauge 
system to the BA equations goes as follows. The inhomogeneities in the BA equation 
correspond to the t'Hooft lines or the singular abelian monopoles. On the other hand 
the Bethe roots correspond to the positions of the BPS monopoles at the particular 
plane. The BA equations reflect the interaction between the t'Hooft lines and the 
BPS monopoles which is S-dual to the process of the W-boson exchange between two 
Wilson lines. 

The second ingredient of the generalized Gaudin concerns the nontrivial symme- 
try breaking in the gauge theory which yields the constant term in the BA equation. 
It is this term that gets renormalized when the Q deformation is added. We expect 
that upon the deformation the theory enjoys the cyclic RG behavior. In the brane 
terms the relevant five- dimensional theory involves the set of M2 branes whose co- 
ordinates are partially fixed by minimizing the superpotential which is generated by 
the string M2 instanton. The positions of the BPS monopoles are fixed by the BA 
equations. 

The relationship between the knot homologies and the cyclic RG seems to be not 
accidental. The point is that cyclic RG can be discovered in the quantum mechanics 
of Calogero -like potential [27] when the Calogero coupling constant is subject of 
renormalization. On the other hand the Calogero model at the rational coupling 
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constant is tired intrinsically to the toric knots. 



6.2 Cyclic RG in Liouville-Toda 



Since there is the AGT relation between the Nekrasov partition function of conformal 
SQCD and the 2d conformal block one could ask about the interpretation of the 
Cooper pairs at the Liouville/Toda side. The relation between the wave functions 
in the Richardson model and the conformal blocks in the perturbed WZW model 
defined on the spectral plane of the fermionic model has been analyzed in [25] . 

Let us briefly comment on the dictionary between the BSC-like model and the 
perturbed conformal theory found in [25J. The Richardson wave function is defined 
by the number of the Bethe roots and can be treated as the conformal block in the 
/3 — 7 system iin the perturbed conformal theory. The coupling constant is introduced 
into the conformal model via the operator 



where <fi(z) is the conventional scalar boson. The Richardson wave function enters 
the integral representation of the perturbed conformal block in the same manner 
as the Gaudin wave function enters the integral representation of the Liouville con- 
formal block (see, for example [39]). The BA equation for the Richardson wave 
function corresponds to the saddle point equation in the integral representation for 
the conformal block. 



The most important observation from the dictionary obtained in [25] is that the 

Cooper pair operator is attributed to the screening operator and the SL(2) algebra 
in the Gaudin model was identified with the algebra of screenings. Semiclassically 
these screening operators are attached to the surface operators hence we have no 
contradiction with the conjectured monopole interpretation. The generalization of 
the dictionary above to the RD model should be similar since the Gaudin model is 
just substituted by the twisted inhomogeneous XXX spin chain. We believe that the 
interpretation of the Cooper pair as screenings works in this well. 

Since the Nekrasov partition function for super conformal QCD is related to the 
conformal blocks in Liouville and Toda theories is is natural to ask what is the pos- 
sible interpretation of the cyclic RG flows in the Liouville/Toda side. We restrict 
ourself by the simple remarks postponing the analysis for the separate study. The 
cyclicity could be looked at in the conformal blocks of the Toda theory corresponding 
to the conformal SQCD. It would correspond to the decoupling of the single flavor 
with the largest mass which enters the conformal dimensions of the vertex operators. 




(6.54) 
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Simultaneously the gauge coupling constant which corresponds to the position of the 
particular vertex operator should be renormalized. One could expect that the Efi- 
mov states in the Liouville/Toda theory could manifest themselves as the particular 
resonant states similar to the example elaborated in [16J. 

6.3 Bion condensates 

If the fermions are identified as the monopoles in the Higgs phase one could con- 
cern on the physical mechanism providing the formation of the bound states with 
monopole charge Qm = 2. The possible mechanism providing the bound states of 
monopoles has been suggested in JUJJ HI] . It is based on the consideration of the 
SYM theory at R 3 x S 1 geometry where the interesting solution with the several 
quantum numbers exists. The most relevant solutions involve both fractional topo- 
logical and magnetic charges. The compact dimension provides the breaking of the 
gauge symmetry and the finite number of the vacuum states. 

In this geometry there are also so-called KK monopoles which provide the close- 
ness of the monopole array. The instanton itself gets interpreted as the bound state of 
the array of monopoles and KK monopole with vanishing magnetic and unit topolog- 
ical charge. The bound state found in |40j involves the monopole with the fractional 
topological charge and the KK-antimonopole with the opposite topological charge. 
Hence the state does not have topological charge at all. 

The bion condensation has been considered in 3+1 dimensions however the bulk 
monopole has the kink counterpart at the string worldsheet. Hence we could con- 
jecture that there are the bound states with the kink charge two similar to the bion 
condensate in the bulk. The naive analysis of the zero modes responsible for the 
attraction supports this possibility however the detailed analysis is required. 

The RG analysis of the model involving the gas of bions and electrically charged 
W-bosons has been considered in |JT] where the RG flows involves the fugacities for 
electric and magnetic components and the coupling constant. The coupled set of the 
RG equations has been solved explicitly in the self-dual case and the solution to the 
RG equations for the fugacities obtained in j?I] is identical to the solution for the 
coupling in the RD model upon the analytic continuation. The period of the RG 
in the solution above is fixed by the RG invariant which has been identified with 
the product of the UV values of the electric and magnetic fugacities y e x y m . The 
similarity between the RG behavior is not accidental since the mapping of the gauge 
theory and the perturbed XY model has been found in [41] . 

In the previous sections we have argued that in the deformed SQCD the period of 
the cycle is fixed by graviphoton field hence one could wonder if the product of fugac- 
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ities has any relation with the graviphoton background. The possible answer could 
be as follows. The presence of the magnetic and electric components in the plasma 
simultaneously implies the possible angular momentum of the field. On the other 
hand it is possible to identify the analogue of the gravimagnetization of the gauge 
theory in the RD model. In the gauge theory it is related to the projection of the 
angular momentum in the Euclidean 4d space-time [45]. The similar differentiation 
of the RD Hamiltonian with respect to e at e = yields 



which is related to the projection of angular momentum as well upon use of the 
commutation relation in SL(2,R). 

Such interpretation suggests the possible place of the full Efimov tower. As 
we have mentioned the new state which appear(disappear) during the RG cycle 
corresponds to the shift in the branch of the logarithm. This is usually attributed 
to the additional flux of the global symmetry charge. Hence the natural conjecture 
is that the higher Efimov states are the bound states of the dyons instead of the 
monopoles. Such states exist only in some region of the moduli space decaying at 
the corresponding curves of the marginal stability. Note that the relation of the 
bion ensembles and the perturbed XY model and interesting interplay of different 
scales has been recently discussed in [22]. We plan to discuss the relation between 
perturbed XY models in [22] and in our paper elsewhere. 

7 Discussion 

In this paper we have focused on the two aspects of the relation between the models 
of truncated BCS-like superconductivity and supersymmetric gauge theories in the 
graviphoton background. It was demonstrated that the vacuum structure of the 
D = 4 bulk theories and D = 2 worldsheet theory at the nonabelian string is 
inherited also by the third dynamical system - some version of the superconducting 
system. The Bethe equations defining the vacuum structure of the gauge theory yield 
the spectrum of the excited states of the Cooper pairs in the superconducting RD 
model. This can be considered as an complicated condition of the self-consistency 
between the several components of the whole system. The consistency between the 
worldvolume theories at D4 and D2 branes with monopole excitation generalizes the 
4d/2d duality for the nonabelian strings. 

The second surprising ingredient of the correspondence concerns the RG behavior. 
It was known for a while that RG behavior in the bulk and worldsheet theories are 




(6.55) 
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consistent. Here we have the third finite-dimensional RD model subsystem which 
enjoys the peculiar cyclic RG behavior. This cyclic behavior can be recognized in 
the D = 2 worldsheet theory and D = 4 bulk theory. In worldsheet theory the 
cycle corresponds to the shift of the branch of the solution. This corresponds to 
the additional flux along the nonabelian string. In the D = 4 theory it can be 
attributed to the change of the branch of the superpotential as well. In the bulk 
theory the change of the branch corresponds to the additional flux string in the 
bulk theory which fits with the transformation of the magnetic string into dyonic 
one. Let us emphasize that the cyclic RG implies the possibility of the multiple 
nonperturbative Aj scales like the multiple gaps in the superconducting model. There 
is some similarity with the scenario with the several nonperturbative scales discussed 
in |3]. 

It is natural to ask if there are other possible manifestations of the cyclic RG 
behavior in the SUSY gauge theories. Let us mention two possible candidates. First 
note that since we have observed the cyclicity in the number of flavors the Seiberg 
duality can be looked at from the new viewpoint. Indeed generically the Seiberg dual 
system has the additional mesonic degrees of freedom and the spectrum is reordered. 
One could assume that the "RG time" in this case is 

t RG = log^f) (7.56) 

and the transition from the initial to the Seiberg dual system could be treated as the 
period of the cyclic RG flow. Naively one needs the second dimensionless parameter 
to get the two-coupling flow. In this case the possible candidate is the ratio of adjoint 
and fundamental masses. Recently the bulk Seiberg duality has been recognized in 
the worldsheet theory [H] where it corresponds to peculiar interchange of the scale 
and orientational moduli of the semilocal nonabelian string. The spectrum in the 
dual worldsheet theory has the additional degrees of freedom as required. 

Secondly it was found recently that the spectrum of the stable BPS particles in 
the worldsheet theory manifests the interesting periodicity [37]. Namely the curves of 
the marginal stability are organized as concentric curves and the transition from one 
curve to the next one corresponds to the shift of the branch of the superpotential. 
The curve of the marginal stability is the very promising place to look at for the 
whole Efimov tower. Indeed as we have discussed above the Efimov states appear 
when the third degree of freedom is added to the two-body system near threshold. 
It is this situation which is realized near the CMS where the bound state of two 
BPS constituents disappears. It seems that phenomena observed in [37] could be the 
manifestation of the RG cycles and we plan to discuss this issue elsewhere. 



20 



The are a few evident directions of the generalization. First, it would be inter- 
esting to generalize the analysis above to the asymmetric XXZ and XYZ twisted 
spin chains corresponding to the 5d and 6d gauge theories [38J. One could expect 
the corresponding superconducting systems for each case where the anisotropy pa- 
rameters should be the particular coupling constants. It would be also interesting to 
discover the superconducting model involving two independent parameters of the Q 
deformation. 

It is evident that the investigation of the cyclic RG flow in the field theories is 
at the very beginning and much more work is required. Among the most immediate 
questions a few could be mentioned. 

• The anomalies can be equally considered as the IR and UV phenomena since 
they have interpretation as the spectrum flow at any scale. This means that 
anomalies are a kind of invariants of the RG cycle. How this could be formu- 
lated in an invariant manner? 

• What is the holographic image of the cyclic RG flow and the cycle step? 

• Is it possible to get the bound state of regulator(UV) and physical (IR) degrees 
of freedom similar to the higher states in the Efimov tower? 

• Is it possible to get the refined partition function involving the Efimov states 
or resonances and is there any relation between such partition functions and 
the knot invariants? The natural candidate to look at is the monodromy of the 
spectrum under the RG cycle. 

We hope to investigate these questions elsewhere. 
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